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Primordial black holes in braneworld cosmologies: astrophysical constraints 
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In two recent papers we explored the modifications to primordial black hole physics when one 
moves to the simplest braneworld model, Randall-Sundrum type II. Both the evaporation law and 
the cosmological evolution of the population can be modified, and additionally accretion of energy 
from the background can be dominant over evaporation at high energies. In this paper we present 
a detailed study of how this impacts upon various astrophysical constraints, analyzing constraints 
from the present density, from the present high-energy photon background radiation, from distortion 
of the microwave background spectrum, and from processes affecting light element abundances both 
during and after nucleosynthesis. Typically, the constraints on the formation rate of primordial 
black holes weaken as compared to the standard cosmology if black hole accretion is unimportant 
at high energies, but can be strengthened in the case of efficient accretion. 



PACS numbers: 98.80.Cq 



I. INTRODUCTION 



The idea that our observable Universe may be a brane 
embedded in a higher-dimensional bulk is one which has 
deep ramifications for cosmology, and which in particu- 
lar may rewrite many of our ideas as to how the Universe 
evolved during its earliest stages. In most cosmological 
contexts, the effects of the braneworld scenario are re- 
stricted to early Universe phenomena, though they may 
impact on present observations by modifying relics from 
the early Universe, such as inflationary perturbations or 
the dark matter density. However, in a recent paper (0, 
hereafter referred to as Paper I), we reported on an excep- 
tion to this rule; Primordial Black Holes (PBHs) formed 
during the early Universe may still probe the bulk di- 
mensions today. 

Our analysis was restricted to the simplest braneworld 
scenario, Randall-Sundrum type II 0, (henceforth RS- 
II), where there is a single bulk dimension of anti-de Sit- 
ter form characterized by an AdS radius of curvature /. 
We showed that provided the AdS radius was sufficiently 
large, PBHs whose lifetime was as long as the present age 
of the Universe could behave as five-dimensional objects, 
and that this would lead to reductions both in the mass 
and the temperature corresponding to a given lifetime. 

In Paper I it was argued that at the time a density 
fluctuation collapses to form a black hole, its mass will 
be of order the horizon mass. Subsequently a braneworld 
black hole could undergo substantial growth by accreting 
material from the cosmological background, as shown by 
Majumdar Q and by Guedens et al hereafter Pa- 
per II). This accretion phase will typically last until the 
standard cosmological regime is reached, after which the 
evolution is governed by Hawking evaporation. 

Once formed, PBHs will influence later cosmological 
epochs, leading to a number of observational constraints 
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on their allowed abundance. These have been extensively 
investigated in the case of the standard cosmology Q- 
[l8[ (see also Ref. for a review) , and the aim of the 
present paper is to reanalyze the main constraints in the 
braneworld context. Because the temperatures of black 
holes evaporating at a given epoch are modified, for the 
most part such constraints have to be recomputed from 
first principles. 

To outline the types of constraints that can arise, let 
us consider an epoch labeled by cosmic time t. PBHs 
whose lifetime exceeds this time will essentially still pos- 
sess their initial mass and only contribute to the overall 
energy density. As the observable Universe is close to flat- 
ness, a conservative bound derives from the fact that the 
PBH mass density should not overdominate the Universe. 
PBHs with lifetimes of order t are evaporating rapidly, 
producing bursts of evaporation products. Limits can be 
obtained from imposing they should not interfere disas- 
trously with established processes such as those of nucle- 
osynthesis, or from the fact that these bursts have not 
been unambiguously observed today 20]. Even shorter- 
lived PBHs will have evaporated completely at an earlier 
stage. If this happened well before the decoupling time of 
a particular species of evaporation product, its Hawking 
radiation will thermalize with the surroundings, boost- 
ing the photon-to-baryon ratio in the process [2l| . In the 
case of evaporation after photon decoupling, the radia- 
tion spectrum remains intact and subsequently redshifts. 
Thus constraints arise from the cosmic ba ckg round ra- 
diation at high frequencies 0, El 0, El H(| ■ Finally, 
if the radiation is emitted in a certain time window be- 
fore photon decoupling, it cannot be fully thermalizcd 
and will distort the CMB spectrum [HlH |24|. If it is 
assumed black holes leave behind a stable relic, this can 
lead to different constraints [25j but will not be pursued 
here. 
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At a given epoch, the impact that is to be constrained 
is usually dominated by those PBHs with a lifetime of 
order the cosmic time at the epoch in question. In the 
accretion phase soon after formation, their energy density 
will fall off more slowly than that of dust. After the accre- 
tion phase we can neglect the energy loss through evapo- 
ration until the last stages of their lifetime. Therefore, in 
the standard cosmological regime the PBHs will predom- 
inantly behave like dust. We conclude that the fraction 
of the total energy due to PBHs grows proportional to 
the scale factor in the radiation-dominated regime, whilst 
staying constant in the matter-dominated regime. Trans- 
lating the observational constraint into an upper limit on 
the initial fraction in PBHs then typically gives extremely 
strong bounds. Furthermore, since the initial BH mass 
is of order the horizon mass, the limit on the initial PBH 
fraction in turn implies a limit on the amplitude of den- 
sity fluctuations, on scales entering the horizon at the 
time of formation of the PBH. 

In this paper we will reconsider constraints from the 
present density of black holes with lifetimes exceeding the 
age of the Universe and from the present photon back- 
ground, as well as the constraint stemming from the lim- 
its on the allowed distortion of the CMB spectrum. We 
will also reconsider constraints arising from the effect of 
PBHs evaporating during or after nucleosynthesis on the 
light element abundances. 



the uncertainty in the non-linear process of black hole 
formation, we introduce a factor / as 



M; = / M H {U) . 



(3) 



In the following sections the constraints will be found to 
be quite insensitive to its precise value. 

The main distinction to be made is whether the PBHs 
are effectively 4D or 5D, which results from comparing 
the event horizon radius ro with the AdS radius I. Inter- 
estingly, with / < 1 a PBH will be small (ro <C /) if and 
only if it formed in the high-energy regime. Similarly, it 
will be large and effectively 4D throughout the bulk of its 
lifetime if and only if it formed in the standard regime. 
The behaviour of such large black holes should reduce to 
that of standard cosmology. 1 

Neglecting possible charges or rotation, the small 
PBHs arc to good approximation described as 5D 
Schwarzschild black holes, for which 
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where an index 4 refers to Planckian quantities as mea- 
sured on the brane. The evolution in time is obtained 
from 



II. THE KEY EXPRESSIONS 
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We begin by reviewing the key results from Papers I 
and II, which can be consulted for the full details. Our 
Universe is taken to be a flat Friedmann brane, with the 
effective 4D cosmological constant set to zero. The en- 
ergy density will be radiation dominated up to the time 
of matter domination in the more recent past. Under 
these conditions, there is an early, high-energy regime in 
which the scale factor, energy density, Hubble radius and 
horizon mass are given in terms of cosmic time t as 



with the accretion and evaporation rates given by 
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Following Paper II, the factor q = AF/n in the ac- 
(1) cretion term parametrizes the efficiency with which the 
black hole accretes the cosmic radiation background, with 
F = corresponding to no accretion and F — 1 corre- 
sponding to perfect accretion of high-energy radiation 



This is followed by a standard regime, in which 
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Here, th is an arbitrary time in the high-energy regime 
and t c = 1/2 is the transition time between the regimes. 

As said, we assume that PBHs form with masses ap- 
proximately equal to the horizon mass. To incorporate 



Recently, an interesting conjecture was made in which the radi- 
ation from RS-II b lack holes can be computed via an AdS/CFT 
equivalence I2(ll27l . If true, this would greatly increase the evap- 
oration rate for large black holes as there are many decay routes 
into conformal field degrees of freedom, potentially leading to 
radically different PBH phenomenology. However, their result 
remains a conjecture, and in this paper we continue to adopt the 
traditional view of PBH evaporation as described in our earlier 
papers QQ. It would be interesting to fully analyze the mod- 
ifications to constraints on PBH abundances using their evapo- 
ration law. 
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from a uniform background. F is therefore expected to 
lie between zero and one. 

In the evaporation rate we have defined 



g « 0.0024 5brane + 0.0012 5bulk , 



(9) 



where gt>ranc is the usual number of particle species, while 
,9buik = 0(1) is the number of bulk degrees of freedom (in 
the simplest case just the five polarization states of the 
graviton) . The energy lost through Hawking evaporation 
is mainly emitted onto the brane. As an example we 
mention the case where the black hole emits only massless 
particles, for which gt rane = 7.25 and g = 0.023. 

Primordial black holes that are relevant for observa- 
tional constraints must have lifetimes greatly exceeding 
the cosmological transition time t c . It is then an excel- 
lent approximation to neglect the evaporation term in 
Eq. until t — t c , and subsequently to neglect the ac- 
cretion term. The mass at the transition reads 
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We stress that this is the mass the PBH has at the ef- 
fective onset of evaporation. The total lifetime t cvap is 
given by the 5D mass-lifetime relation 
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FIG. 1: The initial temperature of black holes evaporating 
at the key epochs of nucleosynthesis (t CVEl p ~ 100s) and the 
present (£ eV a P ~ 14 Gyrs), shown as a function of I. 



I > 10 20 U. As the experimental upper limit on the AdS 
radius currently is quite weak (I < l max « 10 31 Z4 |2qp. 
there is considerable room for 5D PBHs to play a role in 
the cosmological history, including the present. 

A simple reworking of the 5D relations expresses the 
BH mass and temperature at the onset of evaporation in 
terms of I and i e vap : 



By assumption t cvap 3> h,t c . Therefore, for a PBH of a 
given lifetime this relation determines what the mass at 
the onset of evaporation should be, irrespective of the oc- 
currence of accretion. Hence including accretion will not 
change the expected temperature of PBHs evaporating 
at a given epoch. 

The 5D relations above are to be contrasted with the 
usual 4D results, in which accretion plays no significant 
part: 



ro 



Tbh(4D) 
(4D) 



2M 

Ml' 
Mj 
8irM 



1.2 x 10' 



rune 



Mi 



(12) 
(13) 
(14) 



For black holes of a given lifetime i e vap> the question 
arises if they are effectively 4D or 5D. They will be small 
(5D) if the AdS radius is greater than a critical value, 
given by 
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For much smaller values of I, the standard 4D case is 
retrieved. For example, all PBHs with lifetimes up to 
the present age of the Universe would have been five 
dimensional throughout their evolution provided that 
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The temperature of PBHs evaporating today and at 
nucleosynthesis is shown as a function of / in Figure ^ 
Plots of the mass are qualitatively the same. At the 
smallest values of /, the temperature assumes the 4D 
value, for / ks Z m j n (which is equivalent to ro ~ I) accurate 
description is uncertain through lack of exact solutions, 
while for larger I values the BH temperature is reduced. 
Since most of the energy of a PBH is radiated at tem- 
peratures close to the start temperature at the onset of 
evaporation, we conclude that the evaporation products 
present at a certain epoch in cosmology will be cooler 
if I is sufficiently large. For example, PBHs evaporating 
today would produce no massive particles, except in a 
high-energy tail from the late stages of evaporation. 

To conclude, we list the formation times ti in terms 
of a given lifetime t C vap- For a PBH that formed in the 
high-energy regime this is 
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while in the standard regime one finds 
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It is worth bearing in mind that if black holes form 
after inflation, there is an upper limit on the mass scale 
coming from gravitational wave production. In Paper I 
we showed this gives a lower limit on the black hole mass 
at formation of 



Mi > 2 x 10 6 M 5 . 



(20) 



Using the lifetime-formation time relation Eq. (|18p. 
Eq. (|20|l implies a lower limit on the black hole's lifetime: 
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For Z = 10 31 ^4 and 100% efficiency this gives 
t cvap > 1100 fs, 



(21) 



(22) 



i.e. no black holes evaporating until after nucleosynthesis. 
But for an efficiency lower than 86 % or when I/I4 < 10 28 
the lightest black holes allowed by the lower mass limit 
can evaporate before t — Is. Note however that the 
above limit may be a very conservative one. If inflation 
was at an energy scale much lower than the allowed up- 
per limit, the lightest permitted PBHs may evaporate at 
much later times still. 



III. CONSTRAINT FORMALISM 

Constraints on the allowed abundance of PBHs of a 
certain lifetime are formulated as upper bounds on their 
mass fraction. This mass fraction at (Mi) will be defined 
as the ratio of the energy density due to PBHs of initial 
mass Mi and the background radiation density, at a time 
t > ti? 



a t (Mi] 



_ Ppbh,Mi(£) 
Prad(i) 



(23) 



The initial and final mass fractions will be denoted ai 
and a C vap respectively. The purpose of the following sec- 
tions is to reconsider observational constraints on a cva p 
at different cosmological epochs, trace them back to ob- 
tain constraints on the initial mass fraction a- u and to 
compare both types of constraint in the standard and 



braneworld cosmologies. Let us consider a constraint im- 
posed at a given epoch of cosmic time r: 
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in the standard or braneworld scenario respectively. Sev- 
eral constraints imposed at this stage do not depend on 
the individual temperature or mass of the PBHs, but 
only on the total energy contained in them or emitted by 
them, so that Lf T = L™. Examples include the bound 
from the present mass density, from the distortion of the 
CMB, or the deuterium photo-disintegration constraint. 
An example of a bound that does depend on the individ- 
ual characteristics of the PBH is the helium abundance 
constraint (see Section IVII A|) . 

To obtain the corresponding bound on the initial mass 
fraction requires knowledge of the cosmic evolution. Note 
that in the high-energy regime, i.e. in the PBH accretion 
phase, we have 



a t oc M(t) a(t) 



(25) 



whereas in the standard regime the mass fraction simply 
grows with the scale factor (until the final epoch of the 
PBH's lifetime). Tracing Eq. J53J back to the time of 
formation, we get 
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In order to compare the initial constraints Eqs. Ij26(l 
and H27JI . we express their ratio as 
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The ratio of the initial scale factors can be expressed as 
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As mentioned in the introduction, the impact that is to 
be constrained at a certain epoch is usually dominated 
by PBHs with lifetimes of order the cosmic time at that 
epoch. We therefore take r = t cvap . Using Eq. I|10(l for 
M(t c ), and Eqs. {TBI or lO for the lifetime-formation 
time relations, we obtain 3 
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(30) 



2 The reader is warned that different notations are sometimes 
used, such as a = p p bh/ptot- This will not be important in 
the radiation-dominated phase, but the definitions do deviate in 
the matter-dominated phase. 



3 We have omitted powers of /, the ratio of the initial black hole 
mass and the horizon mass (see Eq. 0)- Taking them into ac- 
count we find L™ oc /-(l+2<j)/2(4-<f5 and L ST K /-1/2. Sinc0 
the accretion parameter q is not expected to be larger than 1.5, 
these factors are of peripheral significance. 
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FIG. 2: The behaviour of Eq. il-ilH as a function of the ac- 
cretion parameter q for two choices of the ratio Z/Z m i n , un- 
der the assumption that the imposed limit at evaporation 
is identical in the high energy or standard treatment, i.e. 
icvap = iorap- This shows how constraints are modified com- 
pared to the standard cosmology in situations where the con- 
straint at evaporation is unchanged. If accretion is negligible 
the constraints are weakened in the high-energy case, whereas 
for q > 5/8 accretion leads to the constraints strengthening 
in the braneworld case. 



where Z m i n is the value of the AdS radius at which quan- 
tities obtained in the HE scenario reduce to the stan- 
dard ones, see Eq. I|15fl . This equation is illustrated in 
Fig. [21 Since I > Z m j n , which constraint is the stronger 
is simply determined by the sign of (5 — 8q), at least 
when = L5S„. Thus the initial constraint will be 

stronger in the braneworld scenario if the accretion effi- 
ciency F is more than 49%, whilst becoming weaker for 
efficiencies below 49%. Whatever the accretion efficiency, 
maximum discrepancy with the standard constraint is 



found when / = L 



10°%. Since Z D 



oc t, 



1/3 



, this 



maximum discrepancy rises mildly when considering ear- 
lier epochs, i.e. shorter-lived PBHs. 

Note that if the calculations of Paper II turn out to 
be invalid, and accretion is in fact not an important pro- 
cess even during the high-energy regime, the correct con- 
straints are retrieved by filling in q = in the above 
formalism. 



density exceeds the observed matter density of about 0.3 
of the critical density. Phrased in this way, the constraint 
applies regardless of the presence of a cosmological con- 
stant, and indicates that for any PBHs surviving to the 
present we must have 



a (M) < 



0.3 



7,0 



(31) 



The cosmic microwave background corresponds to a pho- 
ton density as f2 7i o^ 2 = 2.47 x 10~ 5 , with h « 0.7, and 
conservatively we can ignore the cosmic neutrinos. 

From Eqs. (|26|) and (|27|) it is clear that a given obser- 
vational constraint will give the most severe initial con- 
straint for the lightest PBHs it applies to, i.e. those that 
formed the earliest. Employing Eq. I|31|l for PBHs that 
are about to evaporate today, t cvap > to, we find 



* cvap < n 



0.3 



6 x 10 d 
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In the standard cosmology this is a bound on PBHs with 
mass Mi « 2 x 10 19 A/4 = 4 x 10 14 g, and the constraint 
on the initial PBH mass fraction, Eq. (H, reads HEil 
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As examples in the braneworld scenario, we consider 
the case of maximum AdS radius, without accretion 
or with 100% accretion respectively, i.e. the cases I = 
l m „, q = 0, and I = l max , q = 4/-7T, and make use 
of Eq. (J201) ■ The first example gives M { = 10 14 M 4 = 
3 x 10 9 g and 



« < r 1/8 io- 



14 



(34) 



while the second example results in Mi = 6 x 10 M4 
lOg and 



«i < / 



-0.65^q-23 



(35) 



These are examples of the general trend discussed earlier, 
that going to the braneworld case without considering 
accretion weakens the constraint on the formation rate, 
but that including accretion can strengthen it, with fully 
efficient accretion leading to a more powerful constraint 
than in the standard cosmology. 



IV. PBHS MUST NOT OVERDOMINATE THE 
UNIVERSE 

For a particular value of I, Eq. (|18fl or CSt with 
icvap =to = 8 x 10 60 ti gives the formation time of PBHs 
that are evaporating today. Black holes formed later are 
essentially still intact. Their density is constrained by the 
observed matter density in the present Universe; specif- 
ically, we are tracking the relative densities of PBHs to 
radiation, and we must ensure that, given the observed 
radiation density, this ratio does not imply that the PBH 



V. THE PRESENT PHOTON SPECTRUM 

If PBHs evaporate between the time of photon decou- 
pling (idee ~ 10 12 s) and the present day, their radia- 
tion spectra will not be appreciably influenced by the 
background Universe, apart from being redshifted. Thus 
the spectra could constitute a fraction of the cosmic 
background radiation. We will restrict attention to the 
photon spectrum. In Refs. |20l l29l | it was pointed out 
that black holes with temperatures above the QCD scale 
(200 — 300 MeV) should emit quark and gluon jets that 
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subsequently fragment into particles whose rest mass is 
below the black hole temperature. This could signifi- 
cantly alter the spectrum as compared to when the parti- 
cles are emitted directly as Hawking radiation. However, 
for large values of the AdS radius, the temperature of 
PBHs evaporating after decoupling will never be above 
a few MeV, in which case the effect should be negligible. 
We derive an expression for the spectral shape, and cal- 
culate it explicitly in the assumption of a scale invariant 
initial PBH spectrum. We then consider the constraint 
to be imposed on the peak of the spectrum. 



A. The spectral shape 

The spectral photon number emitted onto the brane 
by a small black hole with lifetime t eva , p is obtained from 



l cvap 

dN f a(E) 



dE 



E 2 



2tt 2 exp [E/Tsatt)] - 1 



dt , (36) 



with a(E) the emission cross-section for photons of fre- 
quency oj = E. In terms of the integration variable 



E 



T BH (i) ' 

and using the 5D relations of Section [nj this becomes 



(37) 



dN 



dE 512 7T 4 g \l 4 



M\E- 2 I ' a (x)-£— dx, 
J e x - 1 



( 38 ) 

where Xi = E/T-qh and Tbh denotes the temperature of 
the black hole at the onset of evaporation. In the high- 
frequency limit (E 3> ?bh), all cross-sections reduce to 
the same value (see paper I) 



1 



2^2 



a = 47rr5 = —E l x 

7T 



In this limit, the spectral number becomes 



dN 7T 
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(39) 



(40) 



Note that the spectrum declines as E~ A , compared to the 
E~ 3 tail of the standard spectrum. 

Now consider black holes evaporating at a time t ovap > 
tdcc- We will make the approximation that all the energy 
gets released instantly, but will take the spectrum into 
account. The black hole mass fraction just before evap- 
oration is given by 



M(t c ) a(< evap ) 



(41) 



° vap ~ J M a(ti) 
while the number density is 

/, \ Prad(^evap . ... 

"PBH(Eovap) = Ctcvap jj fu \ ■ \ A1 > 



M{t c 



The energy density in photons of energy E, emitted be- 
tween tevap and £ C vap + ^cvap IS 

dU tov (E) = npBH(tcvap) E 2 ^(tevap) (43) 

U-C/ ^cvap 

We require the present total energy density in Hawking 
photons at a certain energy scale E , denoted as Uq(E q ). 
Radiation emitted at t cvap with energy E will be red- 
shifted to Eq today provided 



E(E ) = E 



a(t ) 

O-(^cvap) 



(44) 



Integrating over all times t evap after decoupling , we have 

to 

U (Eo) = J dU (E ) 

tdcc 

^(^evap) 



^-dU ta AE(E )). (45) 



We will assume ii(to) < t c , i.e. all PBHs light enough 
to have evaporated completely today were formed in the 
high-energy regime. Substituting the relevant formulas 
then results in 4 



ol/4 d \ -4/3 /, N 1/2 
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x / a -. 



MlEl 



, x -59/48 , AT 
tevap \ dN 

U J dE(Eo) 



(^■evap) Q(^evap) 



dt, 



cvap 



cvap 



with the factor Q(t cvap ) incorporating the effect of accre- 
tion. It is defined by 

//\ 27 /4/, N-9/4 

Q( tcva P ) 4(4 - 9)/ " - 2 9 F 9/4 (j-J ( (47) 

and reduces to Q = 1 if accretion is neglected. 

The number spectrum of a black hole of initial tem- 
perature Tbh peaks at an energy E = 6 Tbh, with b w 5 
in the standard treatment |30j . Therefore, unless «i is 
sharply peaked at particular initial epochs, the main con- 
tribution to the integral in Eq. I|46|) is obtained when 
E(Eq) = b Tbh (t cvap), he. from PBHs evaporating at 

^cvap — ^main , where 



E 



b T BH (io) 



12/5 



t . 



(48) 



The contribution from PBHs evaporating earlier will 
come from the high-frequency end of their spectrum, 



4 Having omitted powers of / as discussed in previous sections. 
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while PBHs evaporating at later times will contribute ra- 
diation that originated in the low-frequency end. Using 
the number spectrum Eq. (|3*£)l with Xi = b, we estimate 
the total energy density at energy Eq as 



Uo(Eo) « 0.1 ai .g- 7/1 °^ 69/2 ° t 3 



I \ 1/20 / E 
lH W-u.) -V, 1 . (40) 



_t 4 

29/20 



1/2 /, x -17/10 

to 

U 



with 



4(4-q)/q 



2 9 7r -27/10 6 27/5--18/5 (5Q) 
27/5 / . \ -18/5 / t? \ -27/5 



Eo_ 

Mi 



The shape of the spectrum is seen to depend on the accre- 
tion parameter q. Introducing p through Uq(E ) oc Eq, 
the exponent ranges from p — 1.45 for q = to p = 0.82 
for q = A/tt. For comparison, in the standard scenario 
we have p = 1.5 in this range of the spectrum. 
The time i ma in as defined above will lie in the relevant 
time interval only for energies Eq in the interval 



tdcc 



5/12 



6T BH (to) < E Q < 6T BH (tc 



(51) 



Radiation at lower frequencies will stem completely from 
the low-frequency ends of the instantaneous spectra, with 
the dominant contribution coming from PBHs evaporat- 
ing around idee- Its intensity can generically be neglected 
as compared to the main frequency range. 5 For energies 
Eq > b Tbh(^o)i the dominant part comes from the high- 
frequency tail of PBHs evaporating today. The number 
spectrum Eq. (|38[) with X{ — oo is used to obtain 



Uq(Eq) * W- 3 a ig - 25/16 



1/2 



-41/16 



, x -13/16 , „ x -2 



Mi 



(52) 



B. The observational constraint 

To make contact with observations, it is convenient to 
relate the integrated energy density Uq{Eq) to the spec- 
tral surface brightness I(Eq) through 



I(E ) = 
V °> 4tt Eq 



(53) 



The overall peak in the present spectrum is at £p Ca k 
6Tbh(^o)- Using q = 0, b = 5, t = 8 x 10 6 %, t 



cq 



10 6 to an d 9 = 0.023, we find for I(E pea k), expressed in 
the units keVcm -2 s _1 sr -1 keV -1 : 



-1/16 



(54) 



The strongest constraint to be placed on the PBH spec- 
trum from the observed flux arises at E = Speak- Em- 
ploying the 5D relations of Section [H] for t cvap = to, we 
find that S pea k can range from a hundred MeV to a few 
hundred keV, depending on the AdS radius. Thus if the 
AdS radius is large, PBHs would mainly contribute to 
the hard X-ray background [M El H3 The XRB ex- 
hibits a peak at around 30 KeV, thought to be sourced 
by AGN [34|], although the issue is not settled at present. 
This peak arguably is at too low an energy scale to be ex- 
plained by braneworld PBHs, but note that Teh (to) ~ 50 
keV for I = Z max in our simple set-up. There was long 
thought to be another excess of radiation at 1 to 10 
MeV (the so-called 'MeV bump') [3l],Ei However mea- 
surements made with the COMPTEL telescope [36| have 
shown this to be due to an instrumentation error, and 
the background is now believed to be smooth over a wide 
range of energies. From the steepness of the PBH spec- 
trum Eqs. (|49|l and l|52l) we conclude that, if there is a 
spectrum of radiation due to evaporating PBHs, it should 
be significantly fainter than the observed background. 

As an example we consider I = Z max , to give a peak en- 
ergy £Jp e ak ~ 250 KeV. The observed surface brightness 
at this energy is given by [33| 



*obs 



0.1 keVcm -j! s -1 sr- 1 keV" 



(55) 



The constraint /(-Epcak) < lobs then results in an upper 
limit on the initial mass fraction: 



a; < 10 



-23 



(56) 



For comparison, the corresponding constraint in the stan- 
dard case is obtained from the gamma-ray background at 
Speak ~ lOOMev and reads on < lO" 27 EfllllEllllll- 
An equivalent way to constrain the peak value of the 
photon spectrum is found by expressing the observed 
density of radiation of order a hundred keV in terms of 
its density parameter as 



n « 10" 



(57) 



Assuming that the fraction of the PBH mass going into 
photons is roughly ten percent, one obtains 



fipbh < 10" 



(58) 



which is identical in form to the constraints of Section lTVl 
and simply strengthens them by 8 orders of magnitude. 



VI. DISTORTION OF THE COSMIC 
MICROWAVE BACKGROUND SPECTRUM 



In addition, the Universe is still rather opaque at these early 
times. 



Energy that is released at a time tgz ~ 10 10 to will 
fail to thermalize fully with the background radiation, 
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modifying the Planck law with a chemical potential // 
[^[Si. The injected energy AE is related to /i as 

AE , x 

— = 0.71 /i, (59) 

with E the background energy. Observational results j2^| 
suggest an upper limit on /i given by 

(i < 9 x 1CT 5 . (60) 

This can be used to constrain the fraction of PBHs evap- 
orating around the time tsz, as will now be elaborated. 
Since the period under consideration occurs well after 
neutrino decoupling, only the fraction F of the PBH en- 
ergy released onto the brane that does not go into neutri- 
nos (or gravitons) will be relevant for the above bound. 
For the largest values of the AdS radius, the tempera- 
ture of the black holes evaporating at tsz is such that 
the electron is the only massive particle to be produced, 
leading to an estimate for F as 

Fw0.5. (61) 

Smaller values of the AdS radius presumably increase F 
somewhat, but this will not play a substantial part in our 
order of magnitude estimate of the constraint. When the 
injected energy derives solely from black hole evapora- 
tion, Eq. I|59(l becomes 

Fa cvap {t sz )= 0.71 fi. (62) 

Translating the observational bound Eq. I|60|) into a 
bound on the mass fraction of PBHs evaporating around 
the time tsz results in 

OWap < 1-3 X 10~ 4 = icvap ■ (63) 

For the standard cosmology this gives an initial con- 
straint 

a; < 10 -21 , (64) 

as first obtained in [24| ■ The corresponding initial limits 
in the braneworld scenario for the example of I = l max , 
are E = 10~ 17 for q = and if E = 10~ 28 for q = 4/tt. 

VII. NUCLEOSYNTHESIS CONSTRAINTS 

In principle every cosmological era prior to the time of 
photon decoupling, at around 10 12 seconds, could have 
been affected by PBH particle interactions [lj|. Of these, 
the era of standard big-bang primordial nucleosynthesis 
(SBBNS) is generally regarded as being the best under- 
stood and the most tightly constrained, and so presents 
the best place in which to look for such effects. To 
that end, a number of detailed investigations have al- 
ready provided strong evidence to suggest that a range 
of nucleosynthesis reactions and parameters should in- 
deed have been modified in the presence of evaporat- 
ing PBHs, and that furthermore, for a sufficiently high 



density of PBHs, such modifications would ultimately 
have led to changes in the final light element abundances 
|^l^l^l4Cll4lll42ll43 |. Broadly speaking, in this body 
of work existing observational limits on the light element 
abundances are used to put constraints on the size of 
such modifications, which then typically lead to strong 
constraints on the numbers of PBHs allowed to evapo- 
rate both during and after nucleosynthesis. Here we re- 
examine two nucleosynthesis constraints in the context 
of the braneworld, namely the constraint on the increase 
in production of Hclium-4 due to the injection of PBH 
hadrons j3^, and the constraint on the destruction of 
primordial deuterium by PBH photons |4l| . Our main 
aim here is to get a reasonable estimate of how such con- 
straints are modified, and not to perform a quantitatively 
precise calculation. 

From Eq. i|14|) . one finds that in the standard scenario, 
PBHs which evaporate during nucleosynthesis, which we 
shall take to be between ~ 1 and 400 seconds, have initial 
masses which range from ~ 10 9 g at Is, to ~ 10 10 g at 
400 s, with corresponding temperatures that range from 
- 10 4 GeV to ~ 10 3 GeV. On the other hand, in the 
braneworld picture we find from Eqs. i|16|) and l|17|) that 
PBHs which decay during nucleosynthesis have masses 
at the effective onset of evaporation in the range 

,,(<)-^ s (65) 

with corresponding temperatures in the range 

Strictly speaking, since g is temperature dependent, the 
temperature of the evaporating braneworld PBH has to 
be solved iteratively. However, since the temperature 
dependence of g is in fact very weak over the range of in- 
terest, for most purposes an estimate of g ~ 0.1 usually 
suffices for all temperatures up to a few TeV or so. 6 Ac- 
cordingly, for the maximum value of I allowed by present 
observational limits, I en 10 31 Z4, the PBHs evaporating 
during nucleosynthesis have masses which range from 
~ 10 g at 1 second to ~ 200 g at 400 seconds and have 
corresponding temperatures which range from around 1 
to 0.2 GeV. Thus, we may conclude that at the limit of 
the largest / currently allowed, the 'small' PBHs evapo- 
rating during nucleosynthesis can evaporate a wide range 



6 For standard model fields on the brane and gravity in the bulk 
9brane and 9 respectively range from around 7.25 and 0.023, for 
T < 1 MeV to about 106.75 and 0.26 for 300GeV < T < 1 TeV. 
Therefore, over the temperature range, 1 MeV < T < ITeV, we 
have 0.15 < g 1 / 2 < 0.51 (factor of 4) and 2.9 > ip 1 / 4 > 1.4 
(factor of 2). Hence, it is good enough in most calculations to 
take g ~ 0.1. In principle, however, we note that g could take 
a wider range of values at temperatures much higher than the 
TeV scale. 
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of particles, e.g. massless particles, neutrinos, electrons 
and positrons, muons, pions, etas, kaons, and in addi- 
tion, they are just hot enough to evaporate nucleons. 7 
At the other extreme we can take the smallest value of I 
for which Eqs. I|65|l and (|66|l are valid, i.e. I = Z m i n . Then 
the temperature and mass range of the standard scenario 
arc retrieved, as they must do. 



A. The Helium abundance constraint 

The proton-to-neutron ratio, n/p, is a key parameter 
in primordial nucleosynthesis (for reviews see (4J, l45l lifil 
H l48l I49L l50| ) and small variations in its size can have 
appreciable effects on the final values of the light element 
abundances, 4 He, 3 He, D, 7 Li, that the theory predicts. 

According to the standard picture, n/p 'froze out', 
i.e. fell out of equilibrium and became fixed at a nearly 
constant value, at around 1 second, when the weak in- 
teraction proton-neutron inter-conversion rate, T n ^ p ~ 
G 2 F T 5 , fell below the Hubble expansion rate, H ~ 
(Gat<7*) 1 / 2 T 2 . Since n/p had been kept near Boltzmann 
equilibrium up until this time, at freeze out it would 
have had a value of ~ exp(— AM np /Tf) ~ 1/6, where 
Tf ~ 0.8 MeV was the temperature at freeze-out and 
AM np ~ 1.293MeV the neutron-proton mass difference. 
Subsequently the value of n/p was then affected only by 
neutron beta-decay. In addition, however, neutrons and 
protons at this time were also undergoing collisions and 
were thereby able to form deuterium, via reactions such 
as p + n D + 7. Initially, however, the energy and 

density of the photon background was sufficiently high so 
as to photodissociate all of the D that formed in this way. 
Only after the temperature of the Universe had fallen be- 
low about 0.8 KeV, which occurred at a time of around 
100 seconds, did the photodissociation rate drop below 
the D binding rate and nucleosynthesis start. Once this 
so-called 'D-bottleneck' had been breached, the binding 
of protons and neutrons into D was then quickly followed 
by the subsequent binding of D with protons and neu- 
trons into tritium and 3 He nuclei, and in turn the bind- 
ing of these into 4 He. However, due to the absence of 
both stable mass-five and mass-eight nuclei, and in addi- 
tion the existence of strong Coulomb-barriers to all reac- 
tions that could form nuclei with mass-six, seven, nine or 
heavier, these nucleosynthesis reactions were only able to 
proceed efficiently as far as 4 He. Consequently, the pro- 
cess essentially ended once all of the neutrons (i.e. the 
fuel) had been bound into 4 He nuclei, leaving only a very 
small proportion bound in the form of D and 3 He, to- 
gether with an even smaller proportion that were able to 



overcome the Coulomb barriers and go on further to bind 
into heavier nuclei, mainly 7 Li and 7 Be. 

The production of 3 He and D was so-called rate limited, 
i.e. the precise quantities of these nuclei left over at the 
end of nucleosynthesis were determined by the efficiency 
of their binding into 4 He, or in other words by the reac- 
tion rates of the associated binding processes. Since the 
reaction rates were simply proportional to the values of 
the speed of light, the thermally-averaged cross-sections 
and the density of the baryons, i.e. T a ~ c(o~,T)n a , the 
efficiency of each process would have been sensitive to a 
variation in any of these quantities. Hence, in the absence 
of non-standard physics affecting the cross-sections, the 
determining factor would have been the baryon density, 
rtt,. 8 Thus, the predictions of nucleosynthesis calculations 
are a function of essentially just one parameter, namely 
rib, although this is usually expressed in terms of its pho- 
ton number density normalized form, 77. The greater the 
density of baryons, i.e. the higher the value of 77, the 
faster the reaction rates and the more efficient and com- 
plete the binding of neutrons into 4 He and so the less left 
behind in the form of 3 He and D and vice-versa. 

The production of the heavier nuclei in nucleosynthesis 
are likewise sensitive to 77, but display a more complex 
dependence on its value. 

By the time nucleosynthesis actually started at 100 
seconds, /3-decay had lowered n/p to around 1/7. Given 
that nucleosynthesis ended once virtually all of the neu- 
trons had been bound into 4 He, it therefore follows that 
the mass fraction of He at the end of nucleosynthesis, 
Yp = 4ru He /nb, should have been approximately twice 
the value of the neutron-to-baryon ratio when it began, 
namely 



2n 



2n/p 
+ p 1 + n/p 



0.25. 



(67) 



Although this estimate may seem somewhat crude, it is in 
fact to first order in good agreement with a full numer- 
ical treatment and also with the current observational 
bounds. The production of He has only a mild sensi- 
tivity to 77, since it is guaranteed that virtually all of 
the neutrons will end up in 4 He whatever else happens. 
The small dependence that it does display derives from 
the fact that higher 77 values allow earlier D-bottleneck 
breaching times, implying that nucleosynthesis can start 
sooner and therefore with higher initial values of n/p. 
Aside from the 77 dependence, the actual uncertainty that 
currently arises in the theoretical prediction of Y p , is of 
order 0.2% (cry = 0.0005) and comes primarily from 
the (now small) uncertainty in the neutron lifetime [45j, 
which is presently estimated to be r n = 885.7 ± 0.8s [52j . 



7 With regard to the evaporation of nucleons, it should be also 
noted that, in taking proper account of the grey body factors, it 
is expected that the actual Hawking temperature of these black 
holes will be hotter by a factor of a few than the pure black body 
temperatures quoted here. 



It further follows that, as the baryon density was sensitive to 
the expansion rate of the Universe, any variation in the expan- 
sion rate would also have influenced the final distribution of the 
abundances. 
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During the intervening period between the end of nu- 
cleosynthesis and the present, the universal light element 
abundances are all believed to have undergone a certain 
amount of chemical evolution due to the effects of stellar 
processing. In order to estimate the primordial abun- 
dances today, therefore, it is generally desirable to seek 
out astrophysical sites which have been the least affected 
by this. In particular, since metallicity is generally ex- 
pected to be correlated with the degree of stellar process- 
ing that has taken place |53l IH3 | , sites with low metal- 
licity are thought to be good targets for observation. It 
is then hoped that, with a sufficient understanding of 
the intrinsic physics of such sites, one can extrapolate to 
zero metallicity to obtain the primordial values. The best 
current observational estimates of the cosmic primordial 
4 He abundance are believed to come from studies of he- 
lium and hydrogen recombination lines in low metallicity 
clouds of ionized hydrogen, so-called HII regions, which 
reside in blue compact galaxies |55| . 

Presently, however, estimates of the primordial He 
abundance, extrapolated from observations of these re- 
gions by the different groups of observers, exhibit central 
values that differ by significantly more than their quoted 
statistical errors. This situation is generally believed to 
be indicative of the fact that some or all of these esti- 
mates are dominated by systematics |56j . Though it re- 
mains possible that the scatter could also be evidence for 
a genuine variance in the primordial values themselves. 
In light of this general uncertainty concerning the value 
of Y p , we shall adopt here the compromise value, derived 
from re-analysis of the data by Olive et al. ^3 , who pro- 
posed that 

Yp = 0.238 ±0.002 ± 0.005. (68) 

Here the first error is statistical while the second repre- 
sents an estimate of the overall systematic uncertainty in 
modelling the physics of the HII regions. The errors are 
compatible with all the data reviewed. In the context of 
SBBNS, the above estimate of Y p alone is then consistent 
with 1.2 < T) < 6.3. 9 

As first indicated by Zel'dovich et al. the situation 
just described could have been radically different if one 
allowed for the possibility of a population of PBHs, which 
were hot enough to evaporate nucleons during nucleosyn- 
thesis, i.e. PBHs with Tbh ^ 2GeV, since in such an 
eventuality it turns out (as we shall outline below) that 
n/p would have continued to increase after the weak in- 
teraction freeze-out time, giving rise to the possibility of a 
significantly higher yield of 4 He. Based on observational 
data of the time Zel'dovich et al. calculated that an in- 
crease in n/p of more than about 80% (from its presently 
accepted value), would have resulted in an unacceptable 



Taking account also of the current observational estimates of 
the other primordial abundances yields the smaller 'concordance 
interval' of 2.6 < tjio < 6.2 lil. 



over-production of primordial 4 He. 10 Consequently they 
were able to put a constraint on the mass fraction of 
PBHs evaporating at the time of nucleosynthesis and also 
to translate this into a constraint on the initial mass frac- 
tion. 11 A numerical treatment of the effects suggested by 
Zel'dovich et al., which involved integrating the system 
of nucleosynthesis reaction equations modified by appro- 
priate terms to account for the (annihilation and spalla- 
tion) effects of injected PBH baryons, was later carried 
out by Rothman and Matzner |43. The results of these 
simulations broadly concurred with the semi-analytical 
estimates of Zel'dovich et al. and provided both an im- 
proved quantitative and qualitative understanding of the 
PBH effects they had suggested. 

More recently, however, it has been argued that the 
majority of baryons emitted by PBHs would not be 
emitted directly as nucleons and mesons, as Zel'dovich 
et al. had assumed, but rather via the fragmenta- 
tion of a QCD quark-gluon jet |2^. An analysis of 
this mode of injection, which utilizes the formalism of 
Reno & Seckel |57|. has been carried out by Kohri & 
Yokoyama jl^. Attempting to take proper account of 
QCD effects in this way, they find constraints on the 
mass-fraction in PBHs that are typically one or two or- 
ders of magnitude stronger. 12 An alternative proposition, 
made by Heckler |H,|53, however, is that the QCD jets 
would not directly fragment into hadrons. Instead, the 
particles emitted by a PBH would form a dense plasma, 
through which the quarks and gluons would lose energy 
via QCD bremsstrahlung and pair production, and this 
in turn would give rise to a photosphere, which could 
then be constrained by observations. 

However, since we are here interested in only obtain- 
ing a first semi-analytic estimate of PBH effects in the 
braneworld case, we shall ignore such QCD effects and 
follow the original approach of Zel'dovich ct al. 38J, keep- 
ing in mind its drawbacks and limitations. In order to re- 
consider this constraint for braneworld primordial black 
holes, we shall start by estimating the number of particles 
evaporated by a braneworld PBH during nucleosynthe- 
sis. The total number density of emitted particles, N cm , 
resulting from the complete evaporation of a population 
of PBHs of some given initial mass, may be expressed 



10 Zel'dovich et al. took unacceptable to mean Y p > 0.4. 

11 In addition to studying the effects of PBH injected baryons on the 
value of n/p, Zel'dovich et al. also studied the effects of helium 
spallation by PBH baryons. This they argued would have lead 
to an increase of the deuterium abundance; a claim t hat was 
later supported by the n um erical studies of Vainer et al. I3SI and 
Rothman and Matzner |42| . However, we shall not discuss this 
effect here. 

12 We note that this analysis could in principle be easily extended 
to the present context. Moreover, in contrast to the case of stan- 
dard PBHs, given the somewhat lower Hawking temperature of 
braneworld PBHs evaporating at nucleosynthesis, such an anal- 
ysis would not entail a large extrapolation of the behaviour of 
QCD from presently observed experimental regimes. 
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as 
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iV Bm = 



Ppbh 
{Ecm) ' 



(69) 



where (E om ) is the average energy of the emitted par- 
ticles. The ratio of the energy density in PBHs at 
evaporation to the background radiation energy density, 
Ppbh/Prad, is therefore H3] 



a. 



X< 1 (-Brad) -^rad ' 



(70) 



where (-Brad) and N ra d are similarly the average energy 
and number density of the particles comprising the back- 
ground cosmological radiation fluid. 

To good approximation, the ratio of the average ener- 
gies above is given by the ratio of the PBH temperature 
at the onset of evaporation to the background tempera- 
ture at evaporation. For the four-dimensional PBHs of 
standard cosmology, using Eqs. (fHl and lTH|) and apply- 
ing the standard cosmological temperature-time relation 



(71) 



(see e.g. |4 
time t = t, 
as 381 



the total emitted number density at any 
evap during nucleosynthesis can be expressed 



> X °- 13x .9bra 2 nc3co S 



-1/4 



Mi 

M4 



-1/2 



rad ■ 



(72) 



In the braneworld scenario, using Eqs. (JHJ and (|1 1|) re- 
sults in 



320 



1/4 



M{t c 



-1/2 



|, v ; ,„( — j .'/ .'/,,,, \ ] , j N lad . 

(73) 

Evidently, Eqs. (|73() and 1|72|) have the same functional 
form. However, it is important to recall that both the 
mass and temperature of PBHs of a given lifetime are 
reduced in the braneworld scenario. Taking the mass 
reduction in account, Eqs. I|72(l and (|73|l show that the 
emitted number density for a given mass fraction a CV ap 
is increased in the braneworld scenario. But this is obvi- 
ous from the temperature reduction, which reduces the 
average energy (E cm ) per emitted particle. 

Having found an estimate of the total number density 
of particles evaporated, we now turn our attention to the 
interaction of these particles with the cosmological back- 
ground [3^|. In the simplest case (i.e. in absence of any 



13 Note: In this section we use a notation in which N eln is a number 
density and we reserve n specifically to mean the neutron number 
density. 



PBH chemical potentials) one expects PBHs to emit nu- 
cleons and anti-nucleons with equal measure. First let us 
consider the case of the neutrons and anti-neutrons. An 
anti-neutron emitted by a PBH has two possible fates; 
it may either annihilate with a background neutron, or 
alternatively annihilate with a background proton. At 
high momenta the cross-sections for these processes are 
essentially the same. The first of these possibilities leads 
to no net change in the neutron-to-proton ratio. This fol- 
lows from the fact that on average for every anti-neutron 
emitted by a PBH there is also a neutron emitted, thus 
the background neutron which is annihilated will in ef- 
fect only be replaced on average by another neutron also 
emitted by the PBH. On the other hand, by the same rea- 
soning, in the second interaction the background proton 
is effectively replaced with a PBH neutron. Moreover, 
because there are six times as many protons in the cos- 
mological background as neutrons, this latter reaction is 
six times more likely than the former. Hence, on aver- 
age we expect that for every seven neutrons emitted, one 
simply replaces a background neutron and the other six 
replace background protons, thus increasing the neutron- 
to-proton ratio. 

The story for the protons and anti-protons emitted by 
the PBH is similar. An anti-proton emitted by a PBH 
may either annihilate with a background proton, or with 
a background neutron. Here again it is apparent that the 
first case leads to no net change in the neutron-to-proton 
ratio, as the annihilated background proton will in effect 
just be replaced by a PBH proton, whereas in the second 
interaction, a background neutron is effectively replaced 
with a PBH proton. In this case, however, the latter 
reaction is six times less likely than the former. Hence, 
we expect that for every seven protons emitted, six sim- 
ply replace background protons while only one replaces a 
background neutron. 

To summarize, on average every seven anti-neutrons 
and seven anti-protons emitted by an evaporating PBH 
(along with equal numbers of their anti-particles) effec- 
tively convert twelve background protons and two back- 
ground neutrons into seven protons and seven neutrons. 
Hence, the change in the background neutron number 
density, n Cl is 



A 6 1 

on c = - n em - - p c 



(74) 



where n Gm and p cm are respectively the number densi- 
ties of the neutrons and protons emitted by evaporat- 
ing PBHs, which to a good approximation will be the 
same. Similarly, we find that Sp c w — (5/7)n cm , so that 
S(n+p) = 0. 

As discussed above, current observational estimates 
suggest that the 4 He abundance, Y p should be 23.8 ± 
1.1%, where we have added the errors in quadrature and 
quoted the 2a value. Assuming this to be a legitimate 
conservative estimate, then it obviously follows that the 
largest value that could at present be accommodated by a 
nucleosynthesis scenario that included the effects of PBH 
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evaporations must be 24.9%. The question that we wish 
to answer, however, is how much could PBH effects have 
actually contributed to Y p . To ascertain this we use an in- 
dependent piece of information, namely the lower bound 
on 77 from observations of Lya absorption in quasar spec- 
tra which states that 7710 > 3.4. If we take this 
bound at face value, then since SBBNS conserves 77 and 
PBH evaporations decrease it, it follows therefore that 
nucleosynthesis could not have started at a lower value 
of rj than this. Now for this value of 77 SBBNS predicts 
a value of Y p > 24%. Moreover, since nucleosynthesis 
with PBHs will always produce a higher value of Y p , it is 
impossible to have a nucleosynthesis scenario which in- 
corporates PBH effects for 77 = 3.4 such that Y p will be 
less than this. Thus a necessary condition is that PBHs 
could not have increased the value of Y p by more than 
about 0.9%. Hence, we may estimate that PBHs evapo- 
rating during nucleosynthesis may only increase Y p by as 
much as about 1%, or equivalently that 



2n 



P 



2Sn 



5 2n c 



n + p 7 71 



p < 100' 



(75) 



For all presently allowed values of I, we find that PBHs 
which evaporate during nucleosynthesis emit nucleons. 
In general, however, it is not known exactly what propor- 
tion of the particles injected into the background by the 
PBH will constitute nucleons. In the standard scenario, 
however, Carr has estimated that around 20% of the 
particles emitted by a PBH decaying during nucleosyn- 
thesis will ultimately go into nucleons and anti-nucleons. 
Here, therefore, we shall assume that the PBHs will emit 
a fraction J-N em of the total emitted particles in nucleons 
and anti-nucleons, and assume similarly that T < 0.2. 
Therefore, since n em ~ 3~N eia /4, it follows from Eq. I|75() 
that we must have 



n + p n h 



< 



(76) 



Substituting for N cm using Eq. (j72"|> gives the standard 
constraint 38] 



■tcvap 



< 



0.22 / g, 



1/4 / M- N 1/2 



(77) 



where 77 cvap = n^/N^, is the baryon-to-photon ratio at 
evaporation. If one further assumes that 77 is fixed from 
evaporation to the present day, 1 i.e. 77 ovap = 770, this may 
also be written as 



£Wp < 6 x 10" 



n h h 2 



r \9br. 



1/4 , .,, v 1/2 
.9cos \ I Mi x 



Mi 



(78) 



14 Strictly such an assumption is only applicable in absence of fur- 
ther PBH evaporations after nucleosynthesis, as these would con- 
tribute radiat ion to the background and so lower the baryon-to- 
photon ratio 1611 . However, even if this were to occur, such a 
change is in any case constrained to be fairly small and so should 
not affect our result much. 



where we have used the relation, 770 ~ 2.8 x 10 8 ft\ i h 2 . 

Carrying through the same calculation for braneworld 
PBHs we find 



<W P < 3 x 10 



_ 10 fl h h 2 ( 5co S \ 1/4 ( M(t c 



T V 9 



7i2 



Ma 



1/2 



(79) 



Substituting the factors g — 0.1, #branc = 106.75, g cos = 
10.75, T = 0.2 and n h h 2 « 0.02^ the observational con- 
straints finally are written as |38| 



ttcvap < 1-1 x 10 



,'M-\ 1/2 

-10 I Ju i \ _ ^ST 



(80) 



and 



V Mi 



cvap ' 



Making contact with Section llTTI the ratio of the upper 
limits can be expressed as 



r HE 

cvap 

rST 

cvap 



-1/4 



(82) 



with Z m in/^4 ~ 10 14 — 10 15 for PBHs evaporating during 
nucleosynthesis. By assumption I > Z m i n , hence the con- 
straint will be tighter in the braneworld case. This can 
be understood as follows: The change in helium abun- 
dance due to evaporating PBHs is proportional to the 
emitted number density, see Eq. (|75(l . But for a given 



mass fraction 



this number density is increased for 



braneworld PBHs, as alluded to before. Fig.|21shows how 

the constraints are modified as compared to the standard 
scenario. 

Using the expressions of Section IIIII the bounds im- 
posed at nucleosynthesis can be converted into bounds 
on the initial PBH mass fractions a[. For standard cos- 
mology this gives j33 



an < 3 x 10~ 18 



f M 



1 109, 



-1/2 



(83) 



for PBH masses in the range 10 9 g < M\ < 10 10 g. In the 
braneworld scenario we obtain 



a; < 4.79x10" 



/ ; \ (X-7?)/4(4- 9 ) / x 3(?-l)/4(4-g) 

11 / 1 \ / ^evap \ 



(84) 

with 7j e vap ranging from Is to 400s. 

For the limiting case where I = l max = 10 31 /4, neglect- 
ing accretion (q = 0), and expressing the bound in terms 
of the initial PBH mass Mj this gives 



cti < 3 x 10" 



17 ( M i 

10 g 



-3/8 



(85) 



with Mj in the range lOg < M\ < 200g. 
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FIG. 3: The behaviour of Eq. IHOI as a function of the accre- 
tion parameter q for two choices of the ratio I /Zmin, as appro- 
priate to the Helium constraint. This is to be compared with 
Fig. |21 as the imposed limits at evaporation in the high en- 
ergy and standard treatments, and L^J ap , are no longer 
equal. Their ratio is instead given by Eq. 18211 . In the case of 
the Helium constraint, the limits are more strongly tightened 
in the high-energy regime than for other constraints. Note 
that for the highest values of q, I may not be allowed to be 
quite as large as Z max . 



As remarked in Section [H] for extreme values of the 
AdS radius combined with a highly efficient accretion 
process, PBHs that are to evaporate during the nucle- 
osynthesis era would have been too light to be consis- 
tent with the lower mass limit imposed by inflation. We 
therefore end with an example where accretion is efficient 
(q = 4/7t), but take a more moderate value for the AdS 
radius, I = 10 25 I4. The initial mass fractions are then 
constrained as 



cti < 6 x 10 



-2(i 



M 



10- 



(86) 



for initial masses with range 10 4 g < Mi < 10 2 g. 



B. Deuterium photo-disintegration constraint 

The high-energy particles emitted by evaporating 
PBHs both during and after nucleosynthesis can be suffi- 
ciently energetic to disrupt primordial nuclei. One impor- 
tant reaction of this type is so-called photo- disintegration, 
or in other words, the destruction of primordial nu- 
clei by high-energy PBH photons. Of all the primor- 
dial nuclei deuterium is the most susceptible to photo- 
disintegration, since it has both the highest cross-section 
and also the lowest threshold, Qd ~ 2.25 MeV. A detailed 
analysis of this effect was considered by Lindley |4l| for 
the case of standard four-dimensional PBHs. Here we 
briefly review that work and extend the analysis to the 
context of the RS-II cosmology. 



Evaporation products may interact with the back- 
ground up until recombination. Therefore, in order to 
consider the effects of PBH photo-disintegrations on the 
light element abundances one is generally concerned with 
PBHs evaporating between the end of nucleosynthesis at 
400 seconds and the time of recombination at t Tec w 10 12 
seconds. For standard cosmology, these are PBHs with 
masses in the range 



10 iu g < Mi < lO 13 g 
and temperatures in the range 

10 3 GeV > T BH > 1 GeV. 



(87) 



(88) 



In the braneworld case both ranges will be reduced. For 
when I = l max , the ran 

200 g < M(i c ) < 10 7 g (89) 



the extreme case when I = l max , the ranges become 



and 



200 MeV > T BH > 1 MeV. 



(90) 



Although high-energy PBH photons are directly capa- 
ble of causing photo-disintegration of nuclei, such direct 
disintegrations are in fact extremely rare, since the cross- 
sections for photo-nuclear reactions , even for deuterium 
are extremely small, typically < 10~ 6 . The main ef- 
fect rather is indirect and comes instead from the photo- 
disintegrations caused by the very large number of lower 
energy photons which are produced as a consequence 
of the thermalization of the high-energy PBH photons 
with the background. This thermalization process may 
be understood as follows. The high-energy photons 
emitted by the PBHs interact with the background via 
two main processes, namely via Compton scattering off 
the background electrons and via electron-positron pair- 
production off the nuclei. The energetic electrons and 
positrons produced in these processes then in turn subse- 
quently predominantly lose energy via inverse Compton 
scattering off the background photons. These photons 
then again Compton scatter and pair-produce and so on. 
In this way a single high-energy PBH photon gives rise 
to a 'cascade' of photons, electrons and positrons of in- 
creasing number and decreasing energy. 

The set of photons in each cascade, with energies {Ei}, 
above the threshold Q<\ will form a finite set, and al- 
though these photons will all predominantly interact via 
Compton scattering and pair-production processes, nev- 
ertheless will have a small probability P(Ei) of destroying 
a deuteron. This probability is given by the ratio LJld 
of the mean free path of the deuteron to the photon. If 
one ignores the tiny fractions of 3 He, Li and D, the mean 
free path of a pho ton of energy E in the background is 
approximately |4l| 



Z 7 « n c cr c 



(91) 



where a c and cr pp are respectively the Compton- 
scattering and pair-production cross-sections. Since 
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charge neutrality demands that n p 



this may be ap- 



proximated by I 



at 



deuteron is given by [4l| 



.„pCr t , where the 'total cross-section' 
Hence^ the probability of destroying a 



P(Ei) 



n d a d (Ei) 
n p a t (Ei) ' 



(92) 



where rid and n p are the number densities of the 
deuterons and protons, and a d is the d(-y, n)p cross- 
section. Hence, summing over all the photons in a single 
cascade gives the number of deuterons destroyed, i.e. 



dm 



(93) 



In reality, however, each cascade is different. In order to 
take account of this variance one may introduce an aver- 
age distribution N(E', E), so that the average number of 
deuterons destroyed by a high-energy PBH photon with 
initial energy E is then given by |41| 



dm 



[ N{E',E)P{E')dE'. 



(94) 



In practice the form of N(E',E) can be determined by 
numerically modelling many cascades. Adopting such an 
approach, Lindley was able to empirically estimate 15 the 
above integral to be of the form |41j 



dm 



r>d 







(95) 



Here E* and (3 are constants, with E* ~ 10 1 GeV and 
/3~1. 

Thus far we have focussed only on the effect of photo- 
disintegrations caused by high-energy photons of a single 
energy E. However, PBHs will emit photons with a spec- 
trum of energies. Let us consider a comoving volume V, 
containing a population of PBHs with total mass M(t). 
In a time dt it will evaporate a fraction / 7 dM of its mass 
into photons with a spectrum of energies v(E), such that 



J Ev{E)dE = -f 7 dM. 



(96) 



The number of deuterons in the volume V destroyed in 
this time will therefore be given by integrating over the 
spectrum of cascades that these evaporating photons will 
give rise to, i.e. 



Using Eqs. i|95|l and 11961) we can easily integrate this last 
expression to give |41j 



N d (t 2 ) = X d (t 2 ) 
N d (h) X d (t x ) 



= exp 



AM/ 7 /3 
N P E, 



(98) 



where AM = M(t%) - M{t 2 ) > is the PBH mass evap- 
orated between times t\ and t 2 , N d (ti) is the deuteron 
number contained in V at time ti and X d = n d /rih is the 
deuteron baryon fraction. In principle, Eq. (|98|) should 
strictly apply when PBH evaporation is the only process 
responsible for changing N d , and should therefore apply 
to the time interval between the end of nucleosynthesis 
and recombination, i.e. from ~ 400 seconds until ~ 10 12 
seconds. 16 

Between the times t\ and t 2 deuterons are destroyed, 
so obviously X d (t 2 ) < X d (ti). In addition, in order to 
be consistent with observational limits on the D abun- 
dance, the depletion of X d cannot be too large. Hence, 
we need to demand that X d (t 2 ) > e~ e X d (ti) say, where 
e > is some number of order unity to be constrained by 
observation. 17 With this requirement we therefore have 
the condition 



Q< AMM 
N P E, ~ 



(99) 



Taking iV p M^/m p , where My, is the baryonic mass 
and m p the proton mass, then |41| 



AM e E* 
Mb" ~ 



(100) 



To translate Eq. (|100fl into a bound on some mass frac- 
tion of PBHs, we need to make an assumption about the 
PBH mass spectrum. It is possible PBHs only form in 
a narrow mass range and equivalently only evaporate at 
a specific era. Then t\ and t 2 are simply taken to con- 
tain that era and Eq. (|100fl is effectively a constraint on 
the total PBH mass fraction of the model. On the other 
hand, PBHs could exhibit an extended, smoothly vary- 
ing spectrum. In this case the total mass evaporated is 
dominated by the low mass end of the spectrum, i.e by 
those PBHs evaporating earliest. We may take t\ and t 2 
to be the end of nucleosynthesis and the onset of recombi- 
nation respectively, and estimate AM by the PBH mass 
evaporated shortly after nucleosynthesis. It is therefore 
usually justified to take 



AM 
Mb" 



Ppbh 
Pb . 



(101) 



dN d 



with iovap some time after nucleosynthesis, either when 
N(E' E)P(E') dE' v(E) dE (97) a narrow mass range of PBHs evaporates, or straight 



Lindley also assumed that n^/n c did not change much over the 
thermalization timescale, which is consistent with the underlying 
presumption that nj could not have changed too much. 



In fact the range of validity can be extended back to the begin- 
ning of nucleosynthesis 41], but we shall not consider this here. 
Here Lindley chooses a rough bound of e ~ 1, corresponding a 
decrease in of one e-fold over the period t'2 — ti. 
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after nucleosynthesis for an extended mass spectrum. 
Eq. IjlUUI) is then equivalent to 



-^cvap 



< 



Pb 
Prad 



(102) 



Furthermore, since cc a 3 and p rac i oc a 4 , it follows 
that 



Pb 



Prad 



a(i e 



a(i e 



Pb 



Prad 



o(t) 



t oq a (*eq) 

(103) 

as p ra d = Ptot/2 ~ Pc/2 at the time of matter-radiation 
equality, t cq « 8 x 10 54 i4. The baryon density param- 
eter at equality is readily related to the present one, 
as the matter density parameter at equality is given by 
an(icq) « 1/2: 



ft b ( io ) = ^ m ^ cq ^ b (f u ) w - ^ b ^°^ 
Q m (to) 2 f2 m (to) 

It follows that Eq. 1|1()2J) may also be written as 



*cvap 



< 



m p / 7 /3 



1/2 



»b(*o) 
^m(^o) 



(104) 



(105) 



As a rough estimate for the fraction of the PBH mass 
that decays into photons we take / 7 = 0.1, while the 
ratio of the present baryonic density to the total mat- 
ter density will be taken to be 0.1. A bound on the 
size of the depletion factor, e, can be given by taking 
the difference between the lowest allowed observational 
value and the highest allowed theoretical value of D from 
SBBNS (without PBHs), i.e. assuming the validity of the 
quasar Lyman-aZp/ia bound of rj > 3.4 as before. Follow- 
ing Steigman [g^ , a cautious current observational bound 
on the deuterium abundance is D/H = 3.0lJ'g x 10 -5 . 
On the other hand, a value of r\ > 3.4 implies a maximum 
SBBNS value of D/H < 7.0 x 10~ 5 . Thus, we may take 
e ~ 1, as did Lindley. Filling in all the parameters then 
finally leads to the observational constraint 



< 3.5 x 10 



-29 



1/2 



(106) 



It should be noted that the constraint applies to small 
five-dimensional PBHs produced in the high-energy 
regime as well as to conventional four-dimensional PBHs. 
This is because the amount of deuterium destroyed will 
be proportional to the total amount of emitted energy, 
see Eq. H'J8[1. For a given mass fraction a ova p the emitted 
energy is identical by definition. 

As in the previous sections, we convert Eq. (|1()6|) into 
a constraint on an initial PBH mass fraction, and will 
take tovap = 400s. In standard cosmology this gives 



a>i < 10" 



(107) 



For braneworld black holes, taking I = Z max and q = 0, 
we find 



a s < 10 



-16 



(108) 



Finally, if I = 10 25 /4 and q = 4/tt, the initial constraint 
reads 



-26 



a; < 4 x 10 



VIII. CONCLUSIONS 



(109) 



If there were sufficiently late periods of black hole for- 
mation in the early Universe, their presence could be no- 
ticeable in more recent cosmological epochs such as nucle- 
osynthesis and beyond. Observation puts upper bounds 
on these effects and therefore on the allowed abundance 
of PBHs, conventionally expressed as an upper limit on 
the total PBH energy fraction. Using the evolution equa- 
tions for the black holes and the background cosmology, 
the observational constraints can be translated into con- 
straints on the PBH formation rate. The initial con- 
straints are usually the most severe for those black holes 
with lifetimes comparable with the cosmic time of the 
epoch at which the observational constraint is imposed. 
If PBHs form from the collapse of background density 
perturbations, their formation rate is related to the am- 
plitude of the power spectrum, on scales that enter the 
Hubble horizon around the time the PBHs form. There- 
fore, an initial constraint implies an upper limit to the 
power spectrum on that scale. 

In the RS-II cosmology that we have considered, PBHs 
of a given lifetime i e vap would have formed in the high- 
energy regime, provided the AdS radius of curvature I is 
larger than i m j n oc (tevap/^) 1 ' 3 h- Then the black hole's 
mass and temperature at the onset of evaporation are re- 
duced, and observational constraints have to be adjusted 
in some cases. In addition, and in contrast to four di- 
mensional cosmology, the black hole is likely to grow by 
accretion of the cosmic background as long as it is in the 
high energy phase. Care is needed in this matter, as the 
growth depends very sensitively on the efficiency of ac- 
cretion. Such black holes will be small and effectively 5D 
throughout their lifetime. 

We see that both the PBH and background evolution 
can be altered compared to standard cosmology. As a 
consequence, the translation of an observational limit 
into an upper limit on the initial PBH mass fraction 
will be modified, see Eq. I|30|l . Most of the effects to 
be constrained are simply proportional to the total en- 
ergy in PBHs, rendering the observational constraints in 
standard or braneworld cosmology identical by definition. 
Eq. (|30|) then allows a comparison of the strength of the 
initial constraint in both scenarios. The braneworld con- 
straint will be the weaker if the accretion efficiency is 
below 50%, whilst being stronger for accretion efficien- 
cies above 50%, for all values of the AdS radius I > l m - m - 
It should be noted that in the latter case the initial mass 
of the PBHs is also much smaller than in the standard 
treatment, so that the constraint corresponds to pertur- 
bations on smaller scales. 

If PBHs were sufficiently heavy to have survived to the 
present day, their mass density should not exceed that of 
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dark matter. For PBHs with lifetimes marginally exceed- 
ing to ~ 8 x 10 60 i4, this implies a constraint on the initial 
PBH mass fraction a; . If the PBHs formed in a standard 
cosmological regime, the bound reads ct!j < 10~ 18 . If they 
were formed in the high-energy regime, taking the exam- 
ple of I = Im&x ~ 0.1mm and 100 % accretion efficiency, 
the bounds strengthens to a; < 10~ 23 . 

Primordial black holes evaporating between photon de- 
coupling and the present age leave behind a spectrum 
that peaks at a temperature of order the black hole 
temperature at the onset of evaporation of PBHs with 
iovap ~ ^o- Focussing on the photon component, it is 
required that its density is less than the diffuse cosmic 
background at comparable temperatures. In standard 
cosmology the peak temperature lies around 100 MeV, 
while for I = ^ max it can be as low as 200 keV. In all 
cases, the photon spectrum constraint implies bounds of 
roughly 9 orders of magnitude stronger than the dark 
matter constraint. 

If evaporation products are released around the 
Sunayev-Zel'dovich time tgz — 10 _10 £o> they will fail 
to fully thermalize with the background radiation. This 
time, however, is sufficiently early in order for the excess 
energy to distort the background blackbody spectrum. 
Limits on the allowed distortion of the CMB spectrum 
then imply limits on PBH mass fractions. In standard 
cosmology one obtains ai < 10~ 21 , while the extreme 
braneworld case (I — l max , 100 % accretion) results in 
a; < 10~ 28 . 

If there was a population of PBHs evaporating dur- 
ing or after the era of nucleosynthesis (~ 1 — 400s) , this 
would have led to changes in the final light element abun- 
dances. As a first example, we note that the abundances 
depend on the neutron-to-proton ratio n/p at the on- 
set of nucleosynthesis. A standard calculation predicts 
n/p « 1/7. However, due to PBH evaporation products 
an approximately equal amount of (anti-) protons and 
neutrons is injected into the background, increasing the 
neutron-to-proton ratio. This in turn implies an increase 
in the helium- four mass fraction Y p . Comparing the ob- 
served value with standard theoretical predictions then 
leads to a PBH constraint. Furthermore, the increase 
in Y p is proportional to the number density of evapora- 



tion products. Recall that a PBH of given lifetime has a 
reduced temperature in the braneworld scenario. Thus, 
the associated constraint on the mass fraction a eV a P of 
PBHs at evaporation is strengthened in the braneworld 
case, contrary to what might naively be expected. The 
resulting initial constraint reads a; < 10 -18 in standard 
cosmology. 

For the braneworld case, the following must be borne 
in mind: If there was a period of inflation, an upper 
limit on its energy scale is imposed by the amount of 
produced gravitational waves. The mass of subsequently 
formed PBHs is then bounded from below, in turn im- 
plying a lower limit to the PBH's lifetime. Its strength 
grows with the AdS radius and accretion efficiency, and 
for extreme values PBH evaporation could not have been 
effective until long after nucleosynthesis. We have cho- 
sen I = 10 25 ?4 and 100 % accretion efficiency in the ini- 
tial constraint imposed by helium-four production (Eq. 
(JEU), to obtain a { < 10" 26 . 

In a second example inspired by nucleosynthesis con- 
siderations, it was remarked that evaporation products 
can be sufficiently energetic to destroy newly formed pri- 
mordial nuclei. We have focussed on photo-disintegration 
and considered the change in deuterium abundance, since 
this is the nucleus most susceptible to disintegration. The 
amount of deuterium that would be destroyed grows with 
the total amount of injected PBH energy, showing that 
the ensuing bound on a OV ap does not depend on the in- 
dividual PBH temperature. For the standard initial con- 
straint, one arrives at a; < 10~ 21 . In the example where 
I = 10 25 ?4 and accretion is maximally efficient, this be- 
comes cti < 10~ 26 . 
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